Energy loss in a fluctuating hydrodynamical background 
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Recently it has become apparent that event-by-event fluctuations in the initial state of hydrody- 
namical modelling of ultrarelativistic heavy-ion collisions are crucial in order to understand the full 
centrality dependence of the elliptic flow coefficient 1)2. In particular, in central collisions the density 
fluctuations play a major role in generating the spatial eccentricity in the initial state. This raises 
the question to what degree high Pt physics, in particular leading-parton energy loss, which takes 
place in the background of an evolving medium, is sensitive to the presence of the event-by-event 
density fluctuations in the background. In this work, we report results for the effects of fluctuations 
on the nuclear modification factor Raa in both central and noncentral ^Jsnn = 200 GeV Au+Au 
collisions at RHIC. Two different types of energy-loss models, a radiative and an elastic, are con- 
sidered. In particular, we study the dependence of the results on the assumed spatial size of the 
density fluctuations, and discuss the angular modulation of Raa with respect to the event plane. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The expression "jet tomography" is often used to de- 
scribe the analysis of hard perturbative QCD (pQCD) 
processes taking place inside the soft medium created in 
an ultrarelativistic heavy-ion collision. Jet tomography 
with the aim to study properties of the medium has be- 
come one of the core observables for heavy-ion physics 
at RHIC, and will be even more important in the LHC 
heavy- ion program due to the large kinematic reach of the 
LHC. In particular, the focus so far has usually been on 
the nuclear modification factor Raa between high-energy 
hadron production in A-A collisions and the scaled ex- 
pectation from p-p collisions. The significant suppression 
seen in Raa follows from the energy loss of the leading 
shower-parton caused by its interactions with the soft 
QCD medium, see e.g. (ll-Q- 

Early computations of energy loss were based on rather 
schematic models of the medium, like static cylinders. In 
other words, they implicitly assumed that the mean den- 
sity of the medium is the only relevant tomographical in- 
formation reflected in observables. However, subsequent 
systematic studies of the role of the medium evolution 
model for energy loss 043 have shown that this is not 
the case: Both the azimuth-angle integrated Raa in cen- 
tral collisions as well as the difference between the yields 
of high-px hadrons in the reaction plane and out of the 
reaction plane in non-central collisions reflect the details 
of the medium evolution dynamics [9|. 

On the other hand, in recent years it has been realized 
that details of the bulk-medium evolution, in particular 
the momentum space asymmetries (driven by the dif- 



ferent pressure gradients the in-plane and out of plane 
directions), which are commonly measured as the second 
harmonic coefficient , can only be understood properly 
by taking initial state density fluctuations into account 
[rol - [l2| . This raises a question of how the averaging over 
many events is best performed when the medium is de- 
scribed using relativistic fluid dynamics. The measured 
V2 represents an average over many different events. Pre- 
viously it was implicitly assumed that it is sufficient to 
average over many different initial state geometries and 
run a hydrodynamical simulation once given this aver- 
aged, smooth initial state. However, state of the art 
models now simulate the fluid-dynamical evolution event 
by event (EbyE) for every initial state, and only then 
average over the resulting particle spectra. 

One may therefore also ask to what degree the order 
of averaging matters when considering energy loss in a 
hydrodynamical background medium. An exploratory, 
rather schematic study of this problem has been pre- 
sented in [lH, [l4[ . It is the aim of this paper to study 
the problem in the context of more detailed EbyE fluid- 
dynamical and energy-loss models which are constrained 
by a large body of data. 



II. INITIAL STATE FLUCTUATIONS IN 
HYDRODYNAMICAL MODELING 

We use the event-by-event ideal hydrodynamical 
framework presented in Ref. [Io[ to model the spacetime 
evolution of the bulk QCD-matter produced in ultrarel- 
ativistic heavy ion collisions. 
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A. Initial state 

We apply a Monte Carlo Glauber (MCG) model to 
produce realistic initial states with density fluctuations. 
Nucleons are distributed into each nucleus using a stan- 
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dard two-parameter Woods-Saxon distribution. We have 
not included any finite nucleon-size effects. The collid- 
ing nuclei are separated by an impact parameter b taken 
randomly from a distribution dN / db ~ b. Nucleons i and 
j from the different nuclei arc assumed to collide if 

( ai -« J )» + ( w - Vi )»<™ (1) 

where a^N is the inelastic nucleon-nucleon cross section 
and Xi , yi are the transverse coordinates for the nucleon 
i. 

The initial energy density profile e(x, y) is obtained by 
distributing the density around the transverse positions 
of the wounded nucleons (WN) or binary collisions (BC), 
and using a 2-dimensional Gaussian smearing 

i=l 

(2) 

where a controls the width of the Gaussian. The overall 
normalization constant K, as well as the initial time tq = 
0.17 fm for the hydrodynamical evolution, are motivated 
by the EKRT model [H as in Ref. [10]. For the WN 
(BC) profile K = 37.8 (12.4) GeV/fm. We vary the width 
parameter er in the range 0.4-0.8 fm. 

Centrality classes are defined using the number of WN. 
We slice the distribution of the events into intervals of 
iV-wN so that each interval has a certain percentage of the 
total events. The impact parameter varies freely in each 
centrality class. For simplicity we use the same centrality 
selection for the WN and BC profiles. 

B. Hydrodynamics 

We solve ideal hydrodynamical equations d^T^ = 0, 
where T^ v = (e + P)u^u v ~ g^ v P is the stress-energy ten- 
sor, u 11 is the fluid four-velocity and P is the pressure. 
We simplify the system by assuming longitudinal boost- 
invariance and by neglecting the net-baryon number den- 
sity. Both simplifications are justified since we are inter- 
ested in particle production at mid-rapidity. Pressure is 
related to energy density with an equation of state (EoS). 
Our choice here is the recent lattice EoS from Ref. (l6| . 

Thermal transverse momentum spectra of hadrons 
are calculated with the Cooper-Frye [lTj method from 
a constant-temperature freeze-out hypersurface. We 
choose the freeze-out temperature to be Td ec = 
160(165) MeV for the WN (BC) profile, which leads to 
a reasonable agreement with the measured pion spectra 
in y/SNN = 200 GeV Au+Au collisions at RHIC. In each 
event we have sampled the hadrons from the obtained 
thermal pr spectra, and used PYTHIA [l8[ to do the 
resonance decays as explained in Ref. ficj |. 

When we study the angular dependence of the hard 
partons, we must define the reference plane for each 
event. Theoretically the simplest choice would be the re- 
action plane defined by the impact parameter and beam 



axis. However, in the experiments the impact parameter 
is not known. Thus we follow Ref. [l|| and use the event 
plane as our reference plane 1 . 

The event flow vector Q 2 for each event is defined as 

Q 2 = ^(P Tl cos(2^), iVi sin(2&)), (3) 

i 

where we sum over all particles in the event. The event 
plane is then calculated from the azimuthal angle of the 
vector Q2 as 

arctan(Q y /Q a; ) 

v j ep = ^ ' W 

where arctan is calculated in the correct quadrant. Since 
in our approach we have only a finite number of parti- 
cles in the final state, the event plane fluctuates around 
the " true" event plane, which would be determined from 
infinitely many particles. To eliminate these fluctuations 
we make 200 events from each hydro run and determine 
the final event plane as the average of these 200 event 
planes. This final event plane is now very close to the 
"true" event plane and we can safely neglect the effects 
from the event plane fluctuations. 

The results for final hadron Pt spectra and elliptic flow 
from this EbyE hydrodynamical model with the WN pro- 
file can be found in Ref. [13]. The EoS used in Ref. 
is slightly different, but as shown in Ref. [HI, the differ- 
ences between these EoS with regard to the final state 
hadronic observables are very small. 

For comparison purposes, we also define a "smooth 
initial condition" by averaging over 10 000 MCG con- 
figurations as described above. We do not perform any 
rotation of the configurations used in the averaging pro- 
cess, thus the event plane of the smooth initial condition 
scenario equals to the reaction plane. 



III. ENERGY LOSS MODELLING IN A LUMPY 
BACKGROUND 

When we change from a smooth medium to a medium 
evolving from fluctuating initial conditions, we qualita- 
tively expect to be sensitive to the following four main 
effects on the energy loss of leading partons: 

• In a medium with fluctuations in the initial state, 
the density distribution is not smooth but divides 
into clusters with densities lower and higher than 
the average. Were the suppression of high Pt par- 
tons a phenomenon linear in density, this would not 
matter, but in the observed region Raa is a non- 
linear function of the medium density where the 



Note that the reference plane called "reaction plane" in Ref. figj] 
corresponds to our event plane. 
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effect of reducing the density is stronger than the 
effect of increasing the density by the same factor 2 . 
From these considerations, we can expect the sup- 
pression to weaken when fluctuations (correspond- 
ing to a re-distribution of the smooth medium den- 
sity) are taken into account, i.e. we expect Raa to 
increase. 

• Regions of high and low density in the initial 
state lead early in the hydrodynamical evolution to 
sharp pressure gradients, which immediately start 
to smooth out the fluctuations. This process im- 
plies an irregular flow field in which the local di- 
rection of the velocity field can be quite different 
from the late time radial plus elliptic flow pattern. 
In particular, the flow vector may initially point 
inward to the medium center. To the degree that 
flow influences the strength of energy loss, this ef- 
fect should be visible in the models. A priori, the 
sign of this effect is unknown and needs to be de- 
termined in a calculation. 

• In a fluctuating initial state, clusters of high den- 
sity are distributed event by event around the ver- 
tices of binary collisions (this remains true even if 
wounded nucleon scaling is assumed, as the position 
of a wounded nucleon also implies a binary collision 
at this transverse position). However, since hard 
parton production takes place in binary collisions, 
there is a marked correlation between the initial 
medium density distribution and the hard parton 
production point: Hard partons tend to be pro- 
duced in regions where the matter density is higher 
than average. This increases the suppression and 
is expected to lead to a decrease of Raa when fluc- 
tuations are taken into account. 

• In principle, as discussed in [Tol j . the event plane 
angle differs from the reaction plane angle by a cer- 
tain amount in each event. If the event averaging 
is done relative to the reaction plane rather than 
the event plane, a reduction of the angular mod- 
ulation of Raa may result. Therefore, in the av- 
eraging procedure, we make use of the event plane 
information. 

Note that the hydrodynamical evolution itself 
smoothes the fluctuations over time, i.e. all effects listed 
above become weaker in the later stages of the medium 
evolution. On the other hand, energy loss has a char- 
acteristic length/time dependence inherent to the un- 
derlying physics model. Thus it is not expected that 



2 See e.g. Q] for an explicit computation — a simple argument 
states that Raa is a function bounded from below by zero and 
that complete suppression of partons corresponds to the limit 
of infinite medium density, i.e. a strong increase of density is 
needed to push the suppression below 0.2 



a radiative energy loss model and an elastic energy loss 
model probe the initial state fluctuations in the same 
way. Therefore, we investigate the effect of fluctuations 
in the medium density in two different models for leading- 
parton energy loss. 

In order to illustrate the effect of fluctuations most 
clearly, we use the following strategy: We compute the 
dependence of Raa{Pt = lOGeV, 4>) on the angle </> of 
the outgoing high-P^ hadrons with the event plane in 0- 
10% central collisions for a number N of different events 
with fluctuating initial conditions. For each single event, 
we average over the possibility that the hard parton may 
emerge from any of the binary collision points associated 
with the event. This could be called "ideal tomography" , 
as it contains the maximum possible tomographical infor- 
mation that could be in principle obtained from the event 
(in practice, however, it is astronomically unlikely to find 
a large number of high Pt partons in any single event). 
This results in an R l AA (<f) for the event i which illus- 
trates the intra-event fluctuations as different directions 
lead to paths traversing or passing by dense regions in the 
medium. We compare different R AA (4>) to illustrate the 
magnitude of the inter-event fluctuations, and define the 
final result as the average over the N distinct events as 

Raa{4>) = J2iLo RaaW/N- We determine the medium 
parameters characteristic for the energy loss model by the 
requirement that R AA (^)) / (2n) = i?^ ooth . In other 
words, we require that the result reproduces the angular 
average of a calculation with a smooth initial condition in 
which the medium parameters have been adjusted to de- 
scribe the data at the same Pt- We find no significant Pt 
dependence of the fluctuation effect, which is expected as 
the Pt dependence of Raa is not very sensitive to the 
energy loss model (2fj| . 

We consider three different scenarios for the fluctua- 
tions. In the standard scenario, the initial bulk matter 
distribution follows the WN scaling and the characteris- 
tic size parameter a of the fluctuations is 0.4 fm. In order 
to test the hypothesis that we recover the smooth initial 
condition result for large-sized fluctuations, we study a 
second scenario where a = 0.8 fm. Finally, in order to 
test the sensitivity to the assumption of the WN geome- 
try, we repeat the study by using a BC-distributed energy 
density with a = 0.4 fm. 



A. Radiative energy loss 

We use the Armesto-Salgado- Wiedemann (ASW) for- 
malism in the formulation of [2l| to compute radiative 
energy loss. From a given binary collision vertex and 
with a given orientation with respect to the event plane, 
we compute the two line-integrals along the medium tra- 
jectory 



Q s (r » = (qL) = / (5) 
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and 



B. Elastic energy loss 



uj c (v a ,4>)= / dee $(o 



(6) 



We assume that the transport coefficient q is related 
to the local medium energy density e and the hydrody- 
namically computed local flow velocity p 0] 

q(0 = ifmod • 2 • e 3 / 4 (e)(coshp(e) - sinhp^cosO), 

where a is the angle between the parton trajectory and 
flow vector, and K me d is an adjustable parameter which 
is determined by requiring that the model describes the 
angular-integrated Raa in central collisions. 

Using the numerical results of pl| . Q s and uj c deter- 
mine the energy loss probability density P(AE) which 
we average at a given <f> over all binary collision vertices 
in a given event. As explained in detail in (22[, this pro- 
cedure requires that q(£) along the parton trajectory can 
be written in the form 1 /t q . While this is obviously re- 
alized in the Bjorken model of only longitudinal scaling 
flow, it is not self-evident that this condition is met in a 
hydrodynamical expansion with initial state fluctuations 
and buildup of transverse flow. We have however checked 
for a large number of trial trajectories that the leading 
behaviour is given by l/r a (with a dependent on the 
trajectory and event) and that the deviations are of or- 
der of few percent. This can be understood from the fact 
that the hydrodynamical evolution smoothes out strong 
density fluctuations very early already. 

We denote the resulting average probability density 
as (P(AE)) ( j > . We calculate the momentum spectrum 
of hard partons in leading order (LO) perturbative 
QCD (explicit expressions are given in 0] and references 
therein). The medium-modified perturbative production 
of hadrons at an angle <j) can then be computed from the 
expression 



da 



AA^h+X 



E 



J AA^f+X 

®{P(AE))^D}^ h {z, 4), 

(8) 

where DJ^ h (z, p F ) is the fragmentation function with 
a momentum fraction z at a scale p? F set the hadronic 
Pt [Hj]. From this we compute the nuclear modification 
function R\ A for the event i as a function of the particle's 
angle with respect to the event plane as 



Raa(Pt,V, 



dN l AA /dP T dyd(/> 



(Nbc/<Jnn) daPP/dP T dyd(j) ' 



(9) 



where (Nbc) is the average number of binary collisions 
at a given centrality. 



We model the elastic energy loss of a hard parton as 
discussed in Ref. (24|, by incoherent partonic 2 — > 2 pro- 
cesses in pQCD, with scattering partners sampled from 
the medium. Our simulation of energy losses of high- 
energy partons in the produced QCD matter is based on 
the scattering rate for a high-energy parton of a type i, 



Ti(pi,u(x),T(x)) = 



^2 Fij^. k i(pi,u(x),T(x)), (10) 

j(kl) 



where we account for all possible partonic processes ij —¥ 
kl by summing over all types of collision partners, j = 
u, d, s, u, d, s, g in the initial state, and over all possible 
parton type pairs (kl) in the final state. In general, the 
scattering rate depends on the frame, and in particular 
on the high-energy parton's 4-momentum p\ , on the flow 
4-velocity u(x) and on the temperature T(x) of the fluid 
at each space-time location x. 

In the local rest-frame of the fluid, we can express the 
scattering rate as follows [24| : 



—^J^dE 2f] {E^T)J^ ds[sa, 



(11) 

Here E± is the energy of the high-energy parton i in 
this frame and E% is the energy of the thermal parti- 
cle j with a distribution function fj (E2 , T) , which is 
the Bosc-Einstein distribution for gluons and the Fermi- 
Dirac distribution for quarks. The scattering cross sec- 
tion (Tij_>.fcj( s ) depends on the standard Mandelstam vari- 
able s. A thermal-mass-like overall cut-off scale m = g s T 
is introduced in order to regularize the singularities ap- 
pearing in the cross section when the momentum ex- 
change between partons approaches zero. Here g s is the 
strong coupling constant, which we keep fixed with mo- 
mentum scale in this work. 

To initiate the hard massless parton of a type i in 
each event, we sample the partonic px from the LO 
pQCD single-parton production spectrum da j dprdyi (25j 
at yi — in the range prmin < Pt < \/s/2. For the par- 
ton distribution functions (PDFs), we use the CTEQ6L1 
set [H]. The nuclear effects to the PDFs are con- 

sidered small in comparison with the final state medium 
interactions and thus neglected. The initial rapidity is 
randomly generated in the range \y^\ < j/ max from a fiat 
distribution. This fixes the hard-parton energy E and 
polar angle 6 of its momentum vector p = {px,Py,Pz)- 
The azimuth angle (j>, defined with respect to the event 
plane, is evenly distributed between [0, 2tt]. 

The hard parton is assumed to start interacting with 
the medium at the initial longitudinal proper time to of 
our hydrodynamical model. Since in the cm. frame of 
the colliding nuclei all hard partons are produced in the 
Lorentz-contracted overlap region at z rj 0, the longitu- 
dinal position at later times (before the first collision at 
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t >t ) is assumed to be determined by the longitudinal 
momentum only. The initial time and longitudinal coor- 
dinates for the hard parton are thus to = tq cosh j/j and 
z o = T sinh j/j . The coordinates on the transverse plane 
in the beginning of the simulation are then xq — Xi + %in 
and yo = Vi + if^Q; where the parton position (xi,yi) at 
t = is sampled from the nuclear overlap function for the 
smooth case. With the fluctuating initial conditions, the 
parton starting point is randomly sampled from the set 
of binary collision vertices in each event (see e.g. Fig. [3]). 

The hard parton propagates through the plasma in 
small time steps At, during which we propagate the par- 
ton in position space. The probability for not colliding in 
this time interval is assumed to be given by the Poisson 
distribution e~ riAt , where Ti is the total scattering rate 
(fTTJj) for the hard parton of the type i. For small enough 
At, we can assume that there will be at most one colli- 
sion. As we calculate the scattering rates (fTUf in the local 
rest frame of the quark-gluon plasma fluid element, the 
time step At is also boosted to the same frame. Should 
a scattering happen, the probability Pij^ki for a given 
type of scattering process is determined by the ratios of 
the partial scattering rates (fTT|) to the total scattering 
rate (fit)]) . After scattering, the final state parton with 
highest energy is chosen as the new hard parton to be 
propagated further, for which we repeat the procedure 
outlined above with the next timestep. 

We take into account the system's gradual transfor- 
mation from quark-gluon plasma to hadron gas by us- 

( w \ 1/4 

ing an effective temperature T e ff = ( £ J , where 

9Q = = is the quark-gluon plasma degrees 

of freedom with gluon and quark DOF being g g = 16 
and g q = 6, respectively, and number of quark flavors 
Nf = 3. We always assume there is no significant in- 
teraction between the high-energy parton and the fully 
hadronic medium, and thus no collisions happen in re- 
gions with temperature below Td ec - 

The outcome of the procedure described above is a 
medium-modified distribution of high-energy partons, 
Analogously with the radiative energy loss 



(IN 



AA^f + X 



dp T dy 



case (|8) , the obtained partonic distribution is convoluted 
with the fragmentation function DJ^ h (z, P%) in order 



to 

calculate the nuclear modification factor R aa (Pt , y , </>) ■ 

We average the obtained R aa (Pt , y , (j>) across the ra- 
pidity window [—0.35,0.35] which corresponds to the 
PHENIX acceptance. Due to the non-eikonal propaga- 
tion of the hard partons in the simulation, the initial 
rapidity window is set by choosing y max = 0.7 in order 
to account for all the possible partons falling into the fi- 
nal rapidity window [24] . To achieve roughly the right 
amount of nuclear modification in the 0-10% centrality 
bin and to emulate also the incoherent higher-order pro- 
cesses, the value of strong coupling constant is chosen 
to be a s = 0.5. For the initial conditions in our elastic 
energy loss simulations, the BC profile with the width 
parameter a — 0.4 fm are used exclusively. 



IV. RESULTS 

A. Radiative energy loss, central collisions 

In Fig. [T] we show one of our main results for radiative 
energy loss, i.e. the angular dependence of Raa both for 
a number of events with fluctuating initial state geometry 
in the standard scenario and averaged over 20 different 
events, compared with the result for smooth, averaged 
initial conditions. 
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FIG. 1: The nuclear suppression factor Raa for central 200 
AGeV Au-Au collisions at Pt = 10 GeV as a function of 
the angle of outgoing hadrons with the event plane shown 
with smooth initial conditions, for four different events with 
fluctuating initial conditions averaged over all binary colli- 
sion vertices ('ideal tomography') and averaged over 20 such 
events. 

It is readily apparent from the result that there are 
large inter-event fluctuations as well as somewhat smaller 
intra-event fluctuations. For this centrality class, both 
types of fluctuations are much larger than the angular 
modulation induced by the spatial anisotropy of the mat- 
ter distribution for the smooth initial condition. A 20 
event average captures however the normalization and 
partially also the angular modulation with respect to the 
event plane. For the standard fluctuation scenario, the 
value of -Kmcd in Eq. ((?]) extracted from the data is to 
the limit of our statistical accuracy identical with the 
one from the smooth result, thus we do not observe any 
modification of the medium quenching power in this fluc- 
tuation scenario. 

We may understand this finding better by considering 
the individual mechanisms how fluctuations may influ- 
ence the outcome in detail. In Fig. [2] we show the effect 
of the irregular flow field by comparing the full result with 
a result in which the flow velocity in Eq. (0 has been ar- 
tificially set to zero. This means however that also the 
effect of the mean radial flow field which is present in 
the smooth result has been taken out, thus we rescale 
^mcd in Eq. (JT)) with the factor obtained from a calcula- 
tion with smooth initial conditions where the flow effect 
is also taken out. It is readily apparent from the fig- 
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FIG. 2: The nuclear suppression factor Raa for central 200 
AGeV Au-Au collisions at Pt = 10 GeV as a function of 
the angle of outgoing hadrons with respect to the event plane 
shown for two different events with the flow velocity set to 
zero and evaluating the full flow boost. 



ure that the influence of fluctuation-induced flow on the 
result is small. 
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FIG. 3: The initial state temperature profile, BC points and 
the direction of the event plane for a selected event (num- 
ber 4 in Fig. The constant-temperature contours are for 
T = 160, 300, 400, 500 MeV. Shown as polar plots are also the 
Raa(4>) for this event (solid) and a constant Raa(4>) = 0.22 
(dashed). 



To demonstrate the correlation between the initial 
state density fluctuations and the angular behavior of 
Raa, we have plotted in Fig. [3] the initial temperature 
profile in the transverse plane, the BC points and, as a 



a polar plot, the calculated Raa{4>) for the event num- 
ber 4 of Fig. [TJ We have also added a polar plot of 
Raa{&) — 0.22 to guide the eye. Since many of the 
BC points are in the hot spot in the 1st quadrant, the 
probability to produce a hard parton from this region is 
large. Thus a significant contribution to Raa comes from 
partons emerging from this hot spot. If these partons are 
going into the direction of the event plane, they have less 
medium to traverse and thus do not lose as much energy 
as those which are moving into the opposite direction. 
Hence the Raa is larger near the event plane angle t/'ep- 
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FIG. 4: The nuclear suppression factor Raa for central 200 
AGeV Au-Au collisions at Pt = 10 GeV as a function of 
the angle <f> of outgoing hadrons with respect to the event 
plane shown for two different events averaged with uncorre- 
cted random binary collision points from a smooth nuclear 
overlap function or binary collision vertices correlated with 
the event fluctuations. 



In Fig. 2] we show the angular dependence of Raa av- 
eraged over both randomly chosen initial binary collision 
points according to a smooth overlap distribution and 
averaged over the actual distribution of binary collision 
points of the event. It can be seen from the figure that 
accounting for the correlations of binary collision points 
with medium density leads to increased suppression and 
a reduction of Raa of O(10)% on average. 

These results suggest that the decrease in suppression 
caused by the non-linearity of the response of Raa to the 
medium density is just compensated by the increase in 
suppression due to the correlation of initial vertices with 
dense matter regions, rendering the final result in terms 
of extracted medium properties almost unchanged. How- 
ever, this cancellation is somewhat accidental, as will be 
seen in the discussion of other scenarios for the fluctua- 
tions. 

In Fig. [5] we show the averaged angular modulation 
of Raa for all three fluctuation scenarios discussed pre- 
viously, after a 10 event average where the parameter 
-f^mcd in Eq. ([7]) has been adjusted to result in the best 
fit to the data in the to 10% centrality class, zoomed 
into the relevant y-axis region. Already a 10 event av- 
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FIG. 5: The nuclear suppression factor Raa for 0-10% central 
200 AGeV Au-Au collisions at P T = 10 GeV as a function of 
the angle cf> of outgoing hadrons with respect to the event 
plane shown with smooth initial conditions and as 10 event 
average over fluctuating initial conditions in three different 
scenarios (see text). Note the scale of the y-axis. 



erage recovers some of the angular modulation seen in a 
smooth event. For a — 0.4 fm and WN geometry, the ex- 
tracted .Kmed is, as mentioned above, the same as in the 
smooth case. For a = 0.8 fm we extract a value that is, 
somewhat surprisingly, 18% smaller than in the smooth 
case. Note that this is not a large effect, since the value 
of i^mod may change by a factor two between different 
hydrodynamical evolution models [{J H3| • These results 
suggest that fluctuations do not have a substantial ef- 
fect on the extraction of medium parameters, but rather 
constitute an additional uncertainty. 

The value of if mc d in the BC geometry is 30% larger 
than in the smooth case, but most of this difference can 
be attributed to the density profile geometry rather than 
the fluctuation effect. Vertices for hard parton produc- 
tion are distributed according to the BC geometry. If 
bulk matter is distributed according to the (wider) WN 
geometry, the mean in-medium path for hard partons is 
larger than if bulk matter is also distributed with the 
(narrower) BC geometry. Hence quenching is on average 
stronger for a WN geometry, requiring a larger K mc ^ if 
Raa is required to agree with the data. 



B. Radiative energy loss, non-central collisions 

Let us now turn to the consequences of fluctuations for 
non-central collisions. Here, the angular modulation of 
Raa is much less influenced by fluctuations and already a 
10 event average recovers the sinusoidal modulation with 
a good accuracy. This result is expected and also seen at 
low Pt in the computation of bulk matter V2 ■ 

Our results for the three different scenarios are shown 
in Fig. |U along with PHENIX data for P T = 8.5 GeV and 
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FIG. 6: The nuclear suppression factor Raa for 30-40% pe- 
ripheral 200 AGeV Au-Au collisions at Pt = 10 GeV as a 
function of the angle <f) of outgoing hadrons with respect to 
the event plane shown with smooth initial conditions and as 
10 event averages over fluctuating initial conditions in three 
different scenarios (see text), compared with PHENIX data 

El. 



9.5 GeV pj| . It can be seen that the normalization of 
the results with fluctuating initial conditions is quite dif- 
ferent from the smooth result, although the magnitude of 
the angular modulation is roughly similar. In particular, 
small-sized fluctuations increase the normalization above 
the smooth result, whereas it is reassuring that a larger 
spatial scale for the size of the fluctuations leads back to 
the smooth result. 

Since the Pt dependence of the data is not smooth 
and in particular the 9-10 GeV bin suffers from large 
statistical and systematic errors, it is not entirely clear if 
the data can rule out one scenario at this point. However, 
taking the Py average of the data, it certainly seems 
that the normalization found for small-sized fluctuations 
a = 0.4 fm in the WN scenario is not supported by the 
data. This is an interesting result, as it would allow to 
constrain the typical size of initial state fluctuations. 



C. Elastic energy loss 

Figure [7] is analogous to Fig. [TJ comparing the angular 
dependence of Raa for fluctuating initial state geome- 
try with a result for smooth initial conditions. As the 
individual events do not necessarily correspond to the 
events used in the radiative energy loss model, they are 
here labeled with letters instead of numbers. Also, the 
results are for partons instead of hadrons, as the convolu- 
tion with fragmentation functions requires considerable 
amount of simulation data and increases the statistical 
uncertainties, but contributes very little to the results 
beyond the change in normalization, as can be seen in 
Fig. [7] for the event A. 
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FIG. 7: The partonic nuclear suppression factor Raa for cen- 
tral 200 AGeV Au-Au collisions at pr = 10 GeV as a func- 
tion of the angle of outgoing partons with respect to the event 
plane shown for smooth initial conditions, for four different 
events with fluctuating initial conditions and for an average 
over 20 fluctuation events. Nuclear suppression factor for it 
at Pt = 10 GeV in one event is also displayed. 



Notable inter-event variation was seen in the radia- 
tive energy loss scenario, and this seems to be true also 
with elastic energy loss. The intra-event angular vari- 
ation, however, is rather weak. The average over 20 
events with fluctuating initial conditions, keeping the 
same value a s = 0.5 for the fluctuating and smooth cases, 
equals the smooth initial condition scenario with fairly 
good accuracy, as was the case with the radiative energy 
loss model. This is perhaps more clearly seen in Fig. HI 
in which we also study the scenario where the starting 
points of high-energy partons are sampled from the nu- 
clear overlap function Taa (b) in a medium with fluctuat- 
ing initial conditions, as opposite of using binary collision 
vertices. This effectively removes the local correlation 
between the parton starting points and the regions with 
high energy density in the medium. One would expect 
this to decrease the suppression, and indeed the nuclear 
modification factor increases about 20% in this scenario. 
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FIG. 8: The partonic nuclear suppression factor Raa for cen- 
tral 200 AGeV Au-Au collisions at pr = 10 GeV as a function 
of the angle of outgoing partons with respect to the event 
plane, shown for smooth initial conditions, for an average 
over 20 fluctuation events where parton starting points are 
correlated with binary collision vertices, and for a 20-event 
fluctuation average with parton starting points sampled from 
the corresponding nuclear overlap function TAyt(b). 
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FIG. 9: The partonic nuclear suppression factor Raa for 30- 
40% peripheral 200 AGeV Au-Au collisions at p T = 10 GeV of 
the angle of outgoing partons with respect to the event plane 
shown with smooth initial conditions, and averaged over 20 
events with fluctuating initial conditions. 



Moving to the non-central collisions, we compare Fig.O 
with Fig. [5] As in the radiative energy loss scenario, the 
average of the events with fluctuating initial conditions 
is found to be systematically above the smooth initial 
conditions curve. The difference between the two curves 
is quite small in this scenario, but enough to argue that 
both the radiative and the elastic energy loss models see a 
difference between the smooth and the fluctuating initial 
conditions in the peripheral case. 



V. DISCUSSION 

We have presented a systematic study of the effects of 
QCD-matter density fluctuations on the angular depen- 
dence of the nuclear modfication factor Raa, comparing a 
radiative and an elastic energy loss model. In general, the 
overall effects of fluctuations on Raa observables seem to 
be rather small. 
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In central collisions, an extraction of medium parame- 
ters in the radiative energy-loss model resulted in a dif- 
ference of only ~ 20 % for the different fluctuation 
size-scales studied. For the default fluctuation scenario, 
the 20-event average reproduced both the normalization 
and the angular modulation of Raa of the smooth sce- 
nario with unchanged medium parameter -Kmod- This is 
due a cancellation of the increased suppression due to 
the correlation of parton production points with the hot 
spots, and a decrease of the suppression due the nonlin- 
earity of Raa with respect to the medium density. For 
non-central collisions, however, we find that this cancella- 
tion is incomplete, and as a result the angular-integrated 
Raa depends visibly on the fluctuation size, whereas the 
amplitude of the angular modulation is roughly the same 
as in the smooth scenario. This observation may help 
in constraining the size-scale of the initial QCD-matter 
density fluctuations. 

The results from the radiative and the elastic en- 
ergy loss models seem to agree on qualitative level, even 
though the weak sensitivity of the elastic energy loss 
model to the angle-dependent observables already seen 
in 31] is clearly seen also in this research. In the central 
collisions, no difference is seen between the fluctuating 
and the smooth initial conditions when the average over 
20 events has been taken. However, as in the case of 



radiative energy loss, in the peripheral collisions the fluc- 
tuating conditions appear to produce somewhat smaller 
suppression compared to the smooth background. 

Interestingly, while Fig. [5] would suggest a fairly large 
size for the fluctuations (the a = 0.8 fm case agrees better 
with the PHENIX data), the recent study on the effects 
of fluctuations on thermal photon production [32| would 
favor a smaller fluctuation size. Thus a combined anal- 
ysis has a strong potential to constrain the initial state 
fluctuation scale size and thus to identify the underlying 
physics mechanisms for the bulk QCD-matter produc- 
tion. 
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